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REFRACTION  OF  LIGHT  TRAVERSING  A  CONICAL  SHOCK  WAVE 
PART  I:  CONSTANT  DENSITY,  APPLICATION  TO  SHADOWGRAPHS 


ABSTRACT 

By  means  of  geometrical  optics  the  refraction  of  light  traversing 
a  cone  is  investigated.  It  is  hoped  that  this  study  will  aid  in  the 
interpretation  and  treatment  of  information  from  shadowgraphs  and 
interferograms  used  in  the  investigation  of  the  axially  symmetric 
airflow  about  sharp-nosed  projectiles  in  supersonic  flight. 

In  Part  I,  after  the  necessary  formulas  are  derived,  applications 
to  shadowgraphs  are  made.  (Applications  to  interferometry  will  be 
considered  in  Part  II.) 

A  simple  physical  model  is  assumed:  a  cone  of  index  of  refraction 
located  in  a  space  of  constant  index  of  refraction  n^  (<  Ug)  with 

parallel  light  impinging  on  the  cone  perpendicular  to  its  axis.  The 
main  reason  for  choosing  a  cone  is  that  it  has  simple  geometrical  pro¬ 
perties;  however,  many  of  the  general  ideas  developed  will  be  applicable 
to  any  rotationally  symmetric  Bhock  wave.  In  Parts  I  and  II  a  further 
simplification  is  made;  n^  is  considered  constant.  (Part  III  will  con¬ 
sider  a  variable  appropriate  to  Taylor-Maccoll  flow. ) 

With  Qg  constant  the  ray  tracing  problem  is  fairly  simple  since  no 

differential  equations  must  be  solved.  Snell's  law  and  some  analytic 
geometry  give  the  "exit  functions".  These  designate  the  location  and 
direction  of  light  rays  emerging  from  the  cone  and  are  used  to 
calculate  light  patterns  on  shadowgraphs,  (in  Part  II  they  are  used 
to  calculate  fringe  shifts.)  Since  considerable  labor  is  involved  in 
calculating  them,  the  exit  functions  have  been  computed  on  the  ENIAC 
for  a  wide  range  of  cases,  where  the  parameters  are  the  cone  angle  and 
ratio  of  refractive  indices. 


3 


This  page  Left  Intentionally  Blank 


I  INTRODUCTION 


Optical  methods  are  used  extensively  for  examining  supersonic' flows. 
In  order  to  interpret  the  resulting  photographs  of  the  flow  if  is  necessary 
to  study  the  behavior  of  the  light  which  passes  through  the  supersonic 
disturbance  -  Several  studies  of  the  optical  problems  involved  in  this  in¬ 
terpretation  are  reported  in  references  [lj  .  [2~]  ,  [jj  ,  and  [U]  .  Addi¬ 
tional  references  may  be  found  in  [lQ  .  In  this  report  we  shall  study 
some  of  the  optical  problems  arising  in  the  interpretation  of  shadowgraphs . 

In  Part  II  refraction  effects  in  interferograms  will  be  considered. 

The  schlieren  method  will  not  be  considered  in  detail.  (The  shadowgraph, 
interferogram,  and  schlieren  methods  of  examing  supersonic  flows  can, 
in  their  gross  effects,  be  separated.  However,  there  are  usually  small 
effects  of  one  method  mingled  with  the  gross  effects  of  another.  For 
example,  a  shadowgraph  of  a  free  flight  projectile  was  taken  using  a 
small  rotating  mirror  to  stop  the  motion.  This  shadowgraph  exhibited 
several  fringes  behind  the  shock  wave  and  some  schlieren  effect  where 
one  edge  of  the  mirror  acted  as  a  "knife  edge" . ) 

Consider  the  problem  of  light  ray  tracing  in  the  supersonic  conical 
flow  field  about  the  nose  of  a  cone  cylinder  at  zero  yaw  with  an  attached 
shock  wave.  We  shall  assume  parallel  light  traveling  perpendicular  to  the 
axis  of  the  conical  shock.  The  rays  pass  from  a  medium  of  refractive 
Index  n^,  through  the  conical  region  of  refractive  index  n^(x,y,z) ,  and 

into  the  original  medium  again.  There  they  are  recorded  on  a  photographic 
plate  placed  normal  to  the  initial  beam  of  light,  producing  an  Interfero¬ 
gram  or  a  shadowgraph  according  tc  the  arrangement  of  the  experimental  ap¬ 
paratus  . 

Since  the  change  in  index  of  refraction,  after  the  initial  Jump 
across  the  shock  wave,  is  small  for  large  regions  of  the  flow  field,  it 
is  not  too  unreasonable  to  take  the  index  of  refraction  to  be  constant. 

This  simplifies  the  problem  of  ray  tracing  considerably  since  the  path 
will  be  a  broken  straight  line.  The  difficulty  of  solving  sets  of 
differential  equations  is  eliminated. 

With  the  simplification  of  a  constant  index  of  refraction  the  cal¬ 
culations  required  are  straightforward  but  still  rather  lengthy.  One 
reason  for  this  is  that  the  problem  i3  three-dimensional.  Some  additional 
assumptions  can  be  introduced  to  reduce  the  problem  to  two  dimensions  as 
in  JjsO  .  However,  some  significant  differences  between  the  results  of 
[ 2J  and  the  three-dimensional  treatment  were  found  in  a  preliminary 
investigation.  Hence  the  present  study  was  undertaken. 

From  the  ray  tracing  study  reported  here  qualitative  explanations  can 
be  obtained  of  some  of  the  characteristics  of  shadowgraphs  for  large  values 
of  the  parameters  Mach  number  and  free  stream  density.  The  explanations 
are  qualitative  mainly  because  it  is  necessary  to  consider  rays  other  than 
those  in  a  small  neighborhood  of  the  grazing  ray  so  that  the  constant  in¬ 
dex  assumption  is  not  strictly  justified.  The  calculated  width  of  the 
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shadow  at  a  shock  wave  agrees  fairly  well  (twenty  per  cent)  with  measured 
values.  Equally  good  agreement  can  be  obtained  by  making  the  further 
assumptions  necessary  to  reduce  the  problem  to  a  two-dimensional  one 
(  [2]  ,  [3]  )  except  in  some  cases.  These  cases  are  those  in  which  R/D 
is  either  very  small  or  very  large,  where  R  is  the  radius  of  the  osculating 
circle  of  £2]  and  D  is  the  distance  from  the  disturbance  to  the  photo¬ 
graphic  plate. 


II  DETERMINATION  OF  THE  EXIT  FUNCTIONS 

A  cone  of  refractive  index  n^  is  imbedded  in  a  space  of  refractive 

index  n^,  where  ^4. n^.  The  half -angle  of  the  cone  is  <f.  A  beam  of  rays 

perpendicular  to  the  axis  enters  the  cone,  is  refracted,  and  proceeds  to 
a  photographic  plate.  We  choose  the  coordinate  axes  so  that  z  is  the 
axis  of  the  cone  and  the  rays  travel  in  the  positive  x  direction  (i,j,  and 
k  denoting  unit  vectors  in  the  x,  y,  and  z  directions,  respectively). 

Fig.  2.1  Illustrates  the  coordinate  system  and  the  path  of  the  light. 


Fig.  2.1 


If  r  denotes  the  position  vector  of  a  point  on  the  cone,  the 
equation  of  the  cone  is 

r  .  k  o  z  =  |  r  |  cos  ft , 
or  in  rectangular  coordinates  -- 

x2  +  y2  =  z2  tan2  ft  (2.1) 
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Let  p  be  the  unit  vector  in  the  direction  of  r  ,  i.e.,  p  *  r/  |  r  J  . 
Then  p  •  k  =  cos  (ft. 

The  unit  normal  vector  N  at  the  point  r  lies  in  the  plane 
containing  p  and  k.  From  the  conditions  N  •  N  =  1  and  N  •  p  =  0  the 
outward  unit  normal  is  found  to  be 

N  =  cot  0  p  -  esc  $  k  (2.2) 


At  a  surface  of  discontinuity  in  refractive  index  the  refracted 
ray  lies  in  the  plane  containing  the  incident  ray  and  the  normal  to 
the  surface.  Fig.  2.2  illustrates  the  refraction  of  a  ray  with  angle 
of  incidence  ©1  and  angle  of  refraction  ©2 


n2> 


n. 


(where  is  a  unit  vector  along  the  incident  ray  and  a  unit  vector 

along  the  refracted  ray).  Snell's  law  states  that 

sin02  =  (n^ng)  sin  ^  (2.3) 


Let 


p  =  -  •  N  *  cos 


(2.4) 
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Fig.  2.3  shows  the  projection  of  the  optical  phenomenon  in  the  x, 
y  plane.  We  Bhall  let  z  =  z^-be  the  equation  of  the  incident  plane  beam 

of  light,  and  we  shall  use  p  in  the  equations  as  the  parameter  identi¬ 
fying  the  individual  rays. 


P 


Photographic 

Plate 


(£>Jfc>  >Zp) 


Fig.  2.5 


The  equation  of  an  incident  ray  is  =  1. 


From  eqns.  (2.2)  and  (2.4)  p  =  -  =  -  cot  0  (T  .  p^) 

Since  p1  *  (xj  +  y^3  +  zi^)/(zq  sec  0)> 

P  =  -  (*]_  cot  0)/(z1  sec  0) 

From  this  equation  and  eqn.  (2.1)  we  get  the  coordinates  of  P.  in  terms 
of  p  and  z^  i 


yl 

Z1 


-z^p  tan  0  sec  0 


z1  tan  0  j/l  -  p2  sec2  0 


(2.5) 


Because  of  the  symmetry  we  need  only  study  the  Incident  rays  in  the 
upper  half  of  the  picture,  though  rays  from  the  lower  half  will  be  taken 
into  account  when  we  consider  th  light  pattern  on  the  photographic  plate. 
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As  P1  moves  from  y1  =  z^^  tan  0  (grazing  ray)  to  y^  =  0,  p  variqB  from 
zero  to  cos  0. 


Let  p  be  the  cosine  of  the  angle 
the  emerging  ray  makes  with  the  normal 

i  * 

initial  conditions  p  is  given  by 


(p’)2  =  (1/p2) 


(a  -  Pp)2 


of  exit,  that  iB,  the  angle  that 
at  the  point  Pg.  In  terms  of  the 


/  a2  -  cot2  0  \  2 

\ 1  -  p2  +  cot2  0  / 


(1  -  P2) 


(2.6) 


where 


0  =  n1/n2 

a  =  pp  -  V  p2p2  +  (1  -  p2)'  | 


(2.7) 


For  the  grazing  ray  p  =  0,  and  eqn.  (2.6)  reduces  to 


(p')2  =  -  (Ve2) 


(1  -  P2)2  cot2  0 


(1  -  P2  +  cot2  0)2, 


which  shows  that  p  is  imaginary  and  consequently  the  tangential  ray  exceeds 
the  critical  angle  upon  leaving  the  cone.  One  would  then  expect  a  small 
bundle  of  rays  in  the  neighborhood  of  the  tangential  ray  to  be  totally 
reflected.  (In  all  cases  considered  this  bundle  is  so  small  that  it  has  a 
negligible  effect  on  the  shadowgraph.) 


In  terms  of  initial  conditions  the  coordinates  of  Pg,  the  point  of  exit, 
are  found  to  be  (as  derived  in  the  Appendix) 


where 


x2/zi  =  (l/d)  (tan  0  sec  0)  £ep  cot  0  -  p  (ea  cot  0  +  d)^ 
yg/zj^  =  (l/d)  (tan  0,J  1  -  p2  sec2  0* )  (ea  cot  0  +  d) 


>  (2.8) 


Z2/Z1  =  (Vd) 

d  =  (cot2  0  -  a2)  /  (1  -  p2  +  cot2  0) 
e  -  -  2  (a  -  pp)  cot  0/  (1  -  p2  +  cot2  0) 

The  exit  ray  is  designated  by  the  unit  vector 
=  X  T  +  n  j  +  v  k 


J 


(2.9) 


(2.10) 


# 

See  Appendix  for  derivation. 
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The  direction  cosines  X,  p,  and  v  are 
X 


(l/p)jaf  (e&  cot  0  ‘  p  f®0  cot  0  +  d)}‘  (p°  " 

\i  =  (l/p)  |cos  fijl  -  p2  sec^  w  (eoc  cot  0  +  d)  +  aj 
v  *  (l/p)  [j(a  +  a  )  sin  0j 


J 


where 


*  »  r.  2^  2 

a  =  3 P  -  VU  -  P  )  +  3  P  i 


(2.11) 


(2.12) 


taking  the  positive  root  of  p  in  eqn.  (2.6), 

Eqns.  (2.8)  and(2.1l)  show  that  the  location  and  direction  of  the 
rays  leaving  the  cone  can  be  obtained  by  computing  a  set  of  functions  of 
p  only,  which  will  be  applicable  for  any  incident  beam  of  light  of  the 
form  z  -  constant.  We  shall  denote  (x^/z^),  and  (z 2/zi)  by  Fp  (P)  , 

011,1  fj(p)  respectively  and  refer  to  them  and  X,  n  and  v  as  exit 
functions . 


Exit  functions  were  calculated  on  the  ENIAC  for  a  series  of  cases 

(i.e.,  pairs  of  0  and  p).  The  combinations  of  0  and  B  were  chosen  to 
correspond  to  a  comprehensive  range  of  supersonic  conical  flows. 

Ill  THE  SHADOWGRAPH 


The  image  of  the  shock  wave,  defined  by  the  undeflected  grazing  rays, 
is  its  projection  on  the  photographic  plate  (T^  in  Fig.  3.1,  which  is  a 

cross-sectional  picture).  Consider  a  refracted  ray  which  strikes  the 


Photographic 

Plate 


conical  shock  at  P^.  This  ray  leaves  the  cone  at  P^  with  direction 
cosines  X,  p,  v  and  intersects  a  photographic  plate  located  a  distance  | 


* 

These  functions  are  available  in  the  files  of  the  Airflow  Branch, 
Exterior  Ballistics  Laboratory. 
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(3.1) 


The  functions  G^,  Gg,  G^,  and  were  calculated  on  the  ENIAC  together 
with  the  exit  functions. 

It  is  convenient  to  regard  the  optical  phenomenon  as  a  mapping  of  a 
plane  x  =  constant,  the  object  plane  (Fig.  5.1)#  in  the  uniform  beam  onto 
the  plane  of  the  photographic  plate,  the  image  plane  —  in  other  words,  a 
transformation  from  the  z y^  plane  to  the  z^,  yp  plane. 

Equations  (2 .  ;>)  show  that  the  curve  in  the  object  plane  for  a  constant 

value  of  p  is  y^  =  ^  tan  0^(1  -  p2  sec2  0  V2,  a  straight  line  through  the 

origin,  corresponding  to  a  plane  sheet  of  incident  rays  passing  through  the 
vertex  of  the  conical  shock.  Its  image,  from  eqn.  (3.1)  is 

Vs  =  (W  (zp/fi)  +  (G2G3  -  °iV/°3'  (5‘5) 

a  straight  line  also,  in  general  not  passing  through  the  origin.  In 
particular,  if  the  shock  is  attached  to  a  conical  projectile  of  half  angle 
®  ,  the  locus  of  the  rays  grazing  the  projectile  is  given  by  eqn.  (3.3) 

for  the  value  of  p,  derived  in  the  Appendix,  given  by 

p  =  sec  0  yp2  sin20  -  sin2#  ' /( (3  tan  0)  =  constant  (3.4) 

0  c 

Hence  the  image  of  the  cone  is  a  straight  line  which  does  not 
intersect  the  image  of  the  vertex  of  the  shock  wave,  as  shown  in  Fig.  3.2. 

From  the  computational  standpoint  it  is  convenient  to  study  the 
refraction  of  rays  by  determining  the  images  of  straight  vertical  beams; 
for  moving  from  one  beam  to  another  then  merely  involves  changing  in 

eqn.  (3*1).  Fig.  3.2  presents  the  transforms  of  the  upper  halves  of  a 

series  of  z^  =  const,  curves  for  one  of  the  cases  (0  =  44°,  p  =  .99727; 

O  V* 

with  fl  =  313  )  computed  by  the  ENIAC.  From  this  diagram  one  can  learn 
a  number  of  characteristics  of  a  shadowgraph. 

_ 

This  case  corresponds  roughly  to  a  flow  of  Mach  number  3.2  over  a  cone  of 
half  angle  33>°  at  five  atmospheres  free  stream  pressure  (Fig.  3*4). 
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Note  that  non- dimensional  coordinates  Zp/;|  and  £  are  used. 

These  are  convenient,  since  the  optical  effects  of  placing  the 
photographic  plate  at  different  distances  can  now  he  studied  by 
means  of  one  graph.  Experimentally  the  plate  distance  can  be 
varied  by  using  a  lens  to  obtain  a  " focused  shadowgraph”.  Ordin¬ 
arily,  all  the  phenomena  illustrated  by  a  graph  such  as  Fig.  3.2 
are  not  seen  on  one  shadowgraph.  For  large  §  only  the  part  of  the 
graph  for  small  Zp/fc  is  seen  because  the  large  values  of  Zp/£  are 

actually  out  of  the  conical  flow  region;  for  small  £  the  larger 
values  of  Zp/4  come  into  the  picture,  and  compression  of  the  picture 

causes  loss  of  detail  in  the  region  of  small  Zp /£. 

We  consider  a  particular  curve  of  the  family,  say  z./|  =  .3, 

i  t  i  i  i  1 

and  designate  it  by  K  L  H  G  F  .  (This  is  in  the  range  of  larger 
values  of  z^/S.  For  very  small  z^/5  (<.  .O'j)  the  phenomena  are  somewhat 

different.)  This  curve,  roughly  drawn  in  Fig.  3*3  (b),  is  the  image 

of  EKGF  in  Fig.  3.3(a).  K  corresponds  to  K,  the  ray  with(the  smallest 
value  of  p  that  will  give  a  real  value  (namely,  zero)  to  p  .  As  p 


Vs  */  i 


Fig.  5.3 


increases  (i.e.,  the  object  point  moves  from  K  to  F  In  Fig.  3.3  (a)) 

)  »  «  t  » 

the  image  point  moves  through  K  L  H  G  .  G  corresponds  to  G,  the  ray 
which  grazes  the  conical  body.  All  further  rays,  G  ---  F,  are  blocked 

it 

out  by  the  solid  cone;  thus  the  G  F  segment  of  light  in  Fig.  3.3(b) 

t  f 

does  not  appear  on  the  shadowgraph.  K  L  is  visible,  however;  thus 
some  light  falls  within  the  shadow  of  the  projectile.  This  will  be 
discussed  below. 

The  dark  apace  between  the  projection  of  the  shock  wave  and  the 
z^/ 5  =  const,  curves  is  the  shock  wave  "shadow",  and  its  lower  boundary 
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is  the  envelope  of  these  curves  (MN  in  Fig.  3.2).  It  is  seen  from  the 
behavior  of  these  curves  that  the  inverse  transformation  (yp,zp  to  y^,  z^ ) 

is  double  valued  over  a  large  portion  of  the  region  considered  and  is 
singular  along  this  envelope.  The  shock  wave  shadow  boundary  corresponds 
to  the  locus  of  points  in  the  object  plane  for  which  the  Jacobian  J  = 

d(yp,  zp)/S(yi,  z1)  is  zero.  Setting  the  Jacobian  equal  to  zero  yields 

the  following  relationship  between  z 1  and  p: 

*x(p)  =  5  (g5g2'  -  GjG^)  j  ^a12(o5/o1)'J  , 

where  the  prime  denotes  differentiation  with  respect  to  p.  One  can  then 
in  principle  calculate  the  boundary  of  the  shock  shadow  by  differentiating 
the  G  functions  numerically  and  applying  eqn.  (3.1)  to  give  the  boundary 
curve  in  parametric  form.  However,  in  practice  this  would  be  a  very 
laborious  calculation.  In  Fig.  3.2  the  inner  boundary  of  the  shadow  was 
simply  sketched  in.  Nevertheless,  two  general  conclusions  can  be  drawn 
to  compare  with  the  results  of  [2j  .  (i)  As  Zp/S  approaches  zero  the 

width  of  the  shadow  (measured  perpendicular  to  the  outer  boundary)  does 
not  approach  zero.  Levy’s  approximation  in  [21  gives  a  shadow  width 

l/3 

proportional  to  zp  '  .  (The  image  of  the  tip  region  is  discussed  more 

fully  below.)  (ii)  As  Zp/|  increases,  the  shadow  width  increases  very 

nearly  as  the  one- third  power  of  zp.  But  eventually  this  decreases, 

becomes  zero,  and  thereafter  there  is  no  shadow.  In  [2]  a  monotonic 
increase  is  predicted.  The  fact  that  the  shadow  width  does  not  increase 
monotonically  is  not  illustrated  in  Fig.  3.2.  However,  this  can  be  seen 
by  considering  the  slopes  of  the  straight  lines  p  =  constant  given  by  eqn. 
(3.3).  The  slope  is  always  Gj/G y  and  it  is  found  that  for  the  first  ray 

that  emerges  from  the  cone  (and  a  small  bundle  of  rays  in  the  neighborhood 
of  it)  G^Gj  is  always  greater  than  tan  0.  Hence  eventually  this  bundle 

of  rays  is  imaged  outside  of  the  true  shock  wave.  Thus  for  z p/£  (or  R/D 

in  the  notation  of  [2J  )  small  and  large  the  approximations  of  [2]  are  not 
valid. 

Practically,  to  have  a  large  enough  zp/§  for  a  zero  shadow  width 

it  would  be  necessary  to  have  a  very  small  €  obtained  by  focusing  on  a 
plane  inside  the  disturbance.  In  fact  this  would  be  closely  realized  in 
taking  an  interferogram  where  the  plane  of  focus  is  slightly  off  the 
median  plane.  A  note  of  caution  in  interpreting  Fig.  3*4:  the  zero 
shadow  width  where  the  shock  begins  to  curve  is  not  related  to  the  above 
discussion.  ThiB  arises  because  of  the  expansion  fan  from  the  shoulder, 
and  is  discussed  in  [3J  . 

Further  information  can  be  obtained  from  a  study  of  the  illumination 
of  the  shadowgraph.  Let  1^  be  the  illumination  per  unit  area  (assumed 

constant)  in  the  beam  before  entering  the  disturbance,  and  Ip  (yp,  zp) 
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th  illunlnation  per  unit  area  on  the  photographic  plate.  If  light  from 
the  area  dy^  *  dz^  illuminates  the  area  dyp  •  dzp  on  the  plate,  then 

Ipdypdzp=  dzx  =  ^  [^(yx>  z1)/^(yp,  zp)j  dyp  dzp 

Hence  Ip  =  I1  [  z1)/S(yp,  zp)j  =  I±  J*, 

i 

where  J  is  the  Jacobian  of  the  transformation  yp,  zp  to  y^,  z  .  Actually 

we  should  sum  two  Jacobians  to  account  for  the  fact  that  this  transformation 
is  double  valued.  This  addition  should  take  account  of  the  phase 
difference  also.  However,  we  shall  not  examine  the  interference  effect 
here. 


Since,  as  we  have  Just  seen,  the  Jacobian  J  =  l/j  is  zero  at  the 
boundary  of  the  shock  shadow,  geometrical  optics  predicts  infinite  intensity 
Just  behind  the  shock  wave  on  the  shadowgraph.  Although  this  condition  is 
not  realized  in  actuality,  the  intensity  of  the  illumination  in  this  region 
is  very  high. 

By  locating  the  Images  of  a  series  of  rays  equally  spaced  on  an 
initial  beam  of  light,  one  can  obtain  an  idea  of  the  distribution  of  light 
on  the  shadowgraph.  Thus,  e.g.,  rays  1,  2,  3,  etc.  on  EF  in  Fig.  3*2  are 
taken  at  Intervals  of  l/20  th  of  the  affected  beam  (i.e». ,  l/20  th  of  EF  = 

/  j.,,;  -  ii»  t  i  t  i 

(z1  tan  Their  images  are  points  1,2,3,  etc.  on  K  L  H  F  .  It 

is  seen  that  about  three-fourths  of  the  light  is  blocked  out  by  the  conical 
body.  Only  about  one-fifth  of  the  remaining  light  illuminates  the  segment 

•  i  i 

K  L  H  j  and  since  most  of  this  light  must  be  concentrated  near  the  shock, 

*  t 

most  of  K  L  (rays  for  small  values  of  p)  receives  very  little  light.  Thus 
one  would  still  expect  the  conical  body  to  be  Imaged  sharply,  although  a 
small  amount  of  light  does  enter  its  shadow.  Furthermore,  interference 
fringes  due  to  optical  path  length  differences  between  the  corresponding 
rays  of  two  beams  will  be  noticeable  only  where  both  beams  have  comparable 
intensities,  namely,  near  the  shadow  of  the  shock  wave.  It  can  be  seen 
from  Fig.  3.2  that  for  z^/t  6.0';  light  from  the  upper  half  of  a  beam 

is  imaged  in  the  lower  half  of  the  picture.  Since  most  of  the  image  curve 
corresponds  to  a  very  small_  fraction  of  the  length  of  the  beam,  the 
illumination  in  the  lower  half  of  the  picture  due  to  rays  from  the  upper 
half  of  the  initial  beam  can  4n  large  part  be  neglected. 

The  quantity  1  determines  the  portion  of  the  non-dimensional  picture 
(Fig.  .3-2)  which  is  applicable  to  a  given  shadowgraph.  Thus  for  a  range 
of  very  small  z-J £  a  somewhat  unusual  picture  is  obtained.  The  shadow¬ 
graph  shown  in  Fig.  3.4  -  3;°,  0“  44°,  Mach  no. at  3.2  at  3  atm.) 

involves  conical  flow  only  to  about  z^/§  =  .024,  and  Fig.  3.5  shows  the 
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light  pattern  computed  for  an  approximation  to  this  case.  Only  the  light 
from  the  upper  rays  (y^»0)  which  falls  on  the  screen  is  shown.  The 

small  circles  represent  the  transforms  of  equally  spaced  points  on  the 
z^  =  const,  lines.  The  shock  wave  shadow  by  itself  is  not  visible  here. 

This  shadow  coalesces  with  the  shadow  of  the  conical  body,  and  the  line 
i  r 

0  B  outlined  by  the  light  is  the  image  of  the  conical  body.  Since  a 
very  small  segment  of  the  vertical  light  beam  falls  within  the  shadow, 
and  the  brightest  portion  is  cut  off  by  the  conical  body,  the  illumination 
is  small.  Light  can  also  be  observed  along  0  C  ahead  of  point  C,  the 
point  of  intersection  of  the  cone  shadow  with  the  Zp/§  axis.  The  same 

light  pattern  with  yp  replaced  by  minus  yp  is  produced  on  the  photographic 

plate  by  the  y^<  0  rays. 

Thus  the  image  of  the  boundary  of  the  conical  body  is  outlined  by  light 
rather  than  the  absence  of  light  and  the  nose  or  the  body  is  imaged  into  two 
nappes  of  a  cone.  These  conclusions  are  illustrated  in  Fig.  ;j . >  The  second 
nappe  is  JuBt  barely  visible  because  of  the  small  illumination  as  dis¬ 
cussed  above.  (The  measured  and  the  computed  values  of  OC  (Fig.  5.5) 
differ  by  about  fifteen  per  cent.) 

This  "inversion"  of  the  image  caused  by  the  shadows  coalescing 
and  the  large  deviation  of  rays  near  the  tangential  ray  is  in  general 
verified  by  the  shadowgraph,  Fig.  5>.k.  However,  the  illuminated  region 
inside  the  shadow  of  the  projectile  is  not  of  uniform  intensity;  there  are 
alternate  light  and  dark  bands.  The  following  explanation  is  offered  for 
this  phenomenon:  A  small  roughness  along  the  conical  body  will  produce  a 
weak  curved  shock  wave.  This  will  be  imaged  into  a  dark  band  followed  by 
a  bright  band.  One  side  of  the  image  can  be  calculated  by  a  means  of  eq. 

( 3 • 1 ) *  The  procedure  is  very  similar  to  that  for  finding  the  image  of  the 
conical  body.  One  simply  finds  the  set  of  rays  which  graze  the  wavelet 
(weak  shock)  surface.  This  gives  a  relationship  between. p  and  z^,  and 

the  image  curve  is  then  given  parametrically  by  eq.  (5-l).  In  order  to 
get  the  other  side  of  the  image,  the  rays  which  penetrate  the  wavelet 
surface  would  have  to  be  considered.  It  does  not  seem  necessary  to 
consider  this  further  complication. 

For  simplicity,  the  wavelet  surface  was  taken  to  be  a  cone 
(straight  shock  wave)  forming  a  Mach  angle  with  the  velocity  vector 
at  the  intersection  with  the  main  shock  wave.  This  surface  depends  on 
the  parameters  a  and  where  a  is  the  z  coordinate  of  the  vertex  and 

"  'm~‘  2  2  2  2 

^  is  the  half-angle  (i.e.,  x  +  y  ■  {z-a}  tan  y).  The  set  of  rays 
which  grazes  this  surface  is  given  by 

Va  ■  [-C2  + 
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where 


+  tan2 


solving  a  quadratic 

equation,  the  minus  sign  being  ignored  since  this  would  give  the  rays 
which  graze  the  upstream  nappe  of  the  cone. 

The  image  was  calculated  for  the  conditions  of  case  no.  27  (0  =  44°, 

3  =  . 99727 ),  with  t  =  6l°  33f  and  various  positions  of  the  wavelet.  Th  s 
are  shown  in  Fig.  3.6  for  a/I  =  .01,  .013,  .02,  .03.  As  previously 
remarked,  the  conical  flow  region  extends  only  up  to  about  z/|  =  .024 
in  the  comparison  of  case  no.  27  with  Fig.  3.4.  Thus  only  the  first 
three  values  of  a/|  are  applicable,  and  the  wavelet  images  agree  qualita¬ 
tively  with  Fig.  3-4.  The  dotted  curve,  labeled  "envelope",  is  the 
inner  boundary  of  the  main  shock  wave  shadow,  which  is  very  bright. 

Fig.  3.4  confirms  the  conclusion  that  only  a  small  portion  of  this  bright 
region  is  visible.  Actually,  on  the  basis  of  a  constant  index  of  refraction, 
we  can  calculate  only  one  image,  since  getting  a  second  image  requires 
knowledge  of  rays  that  have  passed  through  a  wavelet.  However,  since  the 
discontinuity  in  index  of  refraction  across  a  wavelet  is  small,  -this  can 
be  neglected  in  calculating  the  images  of  other  wavelets.  Fig.  3*7  is  a 
shadowgraph  taken  under  different  conditions  (M=2.32,  0®33°>  5  atm.  free- 
stream  pressure)  from  those  of  3.4;  it  is  included  because  the  phenomena 
discussed  above  are  more  noticeable. 

In  Fig.  3.6  for  a/ I  =  .03  the  wavelet  image  starts  on  the  image  of 
the  conical  body  and  bends  around  back  into  the  body.  This  is  to  be 
expected  for  the  larger  values  of  z^/l  (smaller  I).  However,  it  would 

be  difficult  to  observe  this  bending  back  because  the  wavelet  image 
would  become  practically  invisible  when  it  is  near  the  bright  region  of 
the  envelope. 

Another  application  of  the  calculations  reported  here  lies  in 
estimating  the  "schlieren  effect"  on  shadowgraphs  and  interferograms 
which  are  taken  by  means  of  a  rotating  mirror  placed  at  the  focus  of 
the  camera  lens.  The  "schlieren  effect"  arises  because  some  of  the 
refracted  light  is  deviated  so  far  from  parallelism  that  it  is  not 
focused  on  the  mirror.  In  an  interferogram  this  causes  the  fringes  to 
b  "washed  out"  in  some  regions.  One  can  use  the  calculated  deviations 
to  determine  the  minimum  size  mirror  which  reduces  the  schlieren  effect 
enough  so  that  good  quality  fringes  are  obtained. 


C1  =  sec2  0  £p2p2tan20  -  32  sin20  sec2^ 
=  tan2  t  jjaPp  tan^0  -  ij 

2  r  2  2  2  2  n 

C5  =  tan  t  j_a  cos  0  sec  ^  +  1  -  a  j 
The  relation  (3.3)  between  z.  and  p  comes  from 
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APPENDIX-DERIVATION  OF  EXIT  FUNCTIONS  AND  SHADOWGRAPH  EQUATIONS 

A  unit  vector  along  the  refracted  ray  inside  the  cone,  1 2,  has  the 
form  (see  Fig.  2.2 ) 

(A-l) 

\  ■  ct  Nx  +  , 

wher  Of  and  0  are  constants  determined  from  eqns.  (2. 5)  and  (2,4),  and 
the  conditions: 

"Ip  •  =  1,  =  -  cos  fig 

The  result  is  p  =  n^/n2  an<^ 

Of  =  Pp  -^pV  +  (1  -  P2)'  (A-2) 

With  eqns.  (2.2),  and  (2.5),  and  (A-l),  and  the  relation  ■  i  one 
can  xpress  in  terms  of  its  cartesian  components: 

*22  -  (p  -  ctp)i  +  (a  cos  fijl  -  p2  sec2  0)j  -  (ot  sin  0)k  (A-5) 

The  refracted  ray  strikes  the  cone  at  P2  (see  Fig.  2.5)  with  position 
vector  r2  and  normal  N^.  Since  r2,  r^,  and  are  coplanar,  p2  ( =  r^/  |  r2  |) 
can  be  written 

Pg  =  +  el2  (A-4) 

where  d  and  e  are  the  constants  defined  in  eqn,  (2.9)  and  determined  from 
the  following  conditions: 

P2  •  P2  =  1,  11  •  k  =  0 

P2  *  k  =  cos  0,  12  •  k  =  ON  ■  k  =  -  a  sin  0 

Since  Jr^J  =  z^  sec  0  and  jr2j  =  z2  sec  0  it  can  be  shown,  using  eqn. 

(A-4),  that  z2/ z1  =  l/d  and 

r2  -  +  (e  zx  sec  0/d)  (A-5) 

From  eqns.  (2.5),  (A-5),  and  (A-5)  we  determine  F^,  F 2,  and  F^.,  the  exit 
functions  listed  in  eqn.  (2.8). 
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If  a  ray  proceeds  from  a  region  of  n 2  to  a  region  of  n^  where  n2 
y  n1,  eqn.  (A-l)  will  give  the  final  ray  here  also.  For  this  simply 
replace  72  by  (-72)  and  by  (-7  ),  and  solve  for  1^. 

\  -  (1/P)  (a*  ^  +  12) 

where  a'  has  the  same  form  as  a,  but  with  a  new  value_of  p,  say  p’.  In 

particular^  consider  the  point  P2  (Fig.  2.3)  and  let  describe  the  ray 

in  n0  and  that  in  n, .  Then 

2  3  1 

7?  =  (l/p)  (a*  N2  +  72)  (A-6) 

where  a*  -  p  p»  -  y/l  -  p2  +  p2  p’2  (A-7) 

and  p’  =7,.  N2;  (A-8) 

i  .e . ,  p*  is  the  cosine  of  the  angle  between  the  emerging  ray  and  the 
normal  to  the  cone  at  the  point  of  exit.  N2  is  obtained  from  eqn.  (2.2). 

N2  =  cot  0  p2  -  esc  0  E  (A-9) 

From  eqns.  (A-6),  (A-7),  (A-8),  and  (A-9)  — 

p7?  .  n2  =  a'  +  72  .  n2 

(p  p'  -  a’ )  *  72  •  n2  *  cot  0  72  .  p2  -  esc  0  72  .  E 
|/fl^  p2)  +  p2  p'2  =  cot  0  72  .  p2  -  esc  0  72  .  E 


From  eqns. 


Since  . 

then 


(A-l),  (A- 3),  and  (A-4)  — 

|/l-  p2  +  P2  p'2  =  cot  0  72  .  (d^  +  e72)  -  esc  0  (-  a  sin  0) 
=  d  cot  0  (72  .  p^)  +  e  cot  0  +  a 
*  (a  +  e  cot  0)  +  d  cot  0  (a  N.^  .  p^  +  .  p1) 

=  (a  +  e  cot  0)  +  pd  cot  0  (7.^  .  f^) 

p^  =  7  .  ^(zi  P  tan  0  sec  0)  T  +  y-jl  +  z-jE^j  /(z-j^  sec  0)  * 

-  -  p  tan  0, 

P2p'2  *  pa  +  e  cot  0)  -  pdpj  2  -  (1-p2) 
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This  last  relation  yields  eqn.  (2.6);  then  a  is  determined  from  eqn. 
(A-7).  N2  can  be  obtained  in  cartesian  form  from  eqn.  (A-9)  ty  applying 

eqns.  (A-4),  (2-5)*  and  (A-5): 


s2  -  r-  pd  +  e  (0  -  ap)  cot  0j  i  +  [(d  +  ea  cot  0)  cos  0 
^0'J  7  +  [d  cot  0  cos  0  -  ea  sin  0  cot  0  - 


•- / ~ 


p  sec 


esc 


0] 


From  eqn.  (A-6)  we  then  obtain  the  direction  cosines  of  the  emerging  ray 
listed  in  eqn.  (2.11). 


The  equation  of  the  line  P2P  in  Fig.  5.1  is 

(  5  -  x2)/x  =  (yp  -  y2)/p  =•  (zp  -  z2)/v, 

where  X,  M,  v  are  given  in  eqn.  (2.11 ).  Remembering  that  F1  <*■  x^/ 

F,  =  y2/Z;L,  F^  =  z2^zl>  one  then  obtains  eqn.  (5.1). 

For  a  shock  wave  attached  to  a  cone  of  half  angle  0 the  first  step 

in  determining  the  image  of  the  cone  is  locating  the  intersection  of  the 
light  rays  with  the  cone.  The  equation  of  the  refracted  ray  within  the 
shock  can  be  expressed,  using 
follows : 

x  =  x±  +  (& 
y  =  yx  +  (c: 
z  =  Zp  ~  (ct 

To  find  the  intersection 

respectively  and  apply  the  equation  of  a  cone,  namely: 

2  ^  2  2t  2, 
x  +  y  =  z  tan  u 
c  J  c  c  c 


eqn.  (A-3),  in  terras  of  the  parameter  s  as 


-  ap)  s 
cos  0V1  -  P  sec  0) 

sin  0)  s 


,  \ 

t)  s  > 


(A-10) 


with  the  cone  set  x,  y,  z  equal  to  x  ,  y  ,  z 


This  leads  to  the  quadratic  equation  in  s 
2 

As  +  Bs  +  C  =  0  r .  -  -  \ 

(A-ll) 

2  2  2 

where  A  =  1  -  a  sin  0  sec  9 

0 

2  2 

B  =  2z1  sin  0  (a  sec  9 -  p0  sec  0) 

C  =  (tan2  0  -  tan20c) 

The  condition  for  the  tangent  ray  is  the  requirement  that  eqn.  (A-ll) 

o 

have  only  one  solution  for  s;  i.e.,  B  -  4AC  =  0.  This  then  yields  eqn. 
(3.4). 
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Diamond  Ordnance  j&tize 
Laboratory 

Connecticut  Avenue  and 

Van  Ness  Street,  N.  W.  1 
Washington  ?.  j}  D.  C. 

Of  Interest  to:  Dr.  Irvine 
C .  Gardner 

*'■  Dr  .  L .  Marton 


Applied  Physics  Laboratory  , 

Tho'ctl  Georgia  Avenue  h*  WP**)  / 

Silver  Spring,  Maryland 


American  Optical  Company 
14  Mechanic  Street 
Southbridge,  Massachusetts 


Farrand  Optical  Company 
Bronx  Blvd.  and  2J-8  Street 
New  York  70,  New  York 

Massachusetts  Institute  of 
Technology 

Department  Mechanical 
Engineering 

Cambridge  ;>9,  Massachusetts 
Attn:  Professor  A.  H. 

Shapiro 

James  Forrestal  Research 
Center 

Princeton  University 
Princeton,  New  Jersey 
Of  Interest  to:  Prof. 

S.  Bogdonoff 


Bausch  and  Lomb  Optical  1 

Company 

262  St.  Paul  Street 
Rochester  2,  New1  York 

Brown  University  ± 

Division  of  Applied  Mathematics 
Providence  12,  Rhode  Island 
Attn:  Professor  R.  F.  Probstein 

California  Institute  of  1 

Technology 

Guggenheim  Aeronautical 
Laboratory 
Pasadena,  California 
Attn:  Professor  L.  Lees  ^ 

Cornell  Aeronautical  Laboratory 
Genesee  Street 
Buffalo  New  York 
Attn:  Miss  Elma  Evans  ^ 

Eastman  Kodak  Research 
Laboratory 

Rochester  ly,  New  York  ^ 

Attn:  Mr.  M.  Herzberger 
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r  Massachusetts  Institute  of 
Technology 

Gas  Turbine  Laboratory 
Cambridge  29,  Massachusetts 
Attn:  Prof.  Ernest  Neumann 

.  ^  North  American  Aviation, 

Inc . 

Aerophysico  Laboratory 
Los  Angeles  k^,  California 

Lehigh  University 
Department  of  Physics 
Bethlehem,  Pennsylvania 
Attn:  Prof.  R.  J.  Emrich 

Perkin- Elmer  Corporation 
y.iy  Hope  Street 
Glenbrook,  Connecticut 
Attn:  Mr.  J.  Baker 

Polaroid  Corporation 
Cambridge  ;9,  Massachusetts 
Attn : Mr .  David  Grey 

Purdue  University 
Department  of  Mechanical 
Engineering 
Lafayette ,  Indiana 
Attn:  Mr.  Vollmer  E. 
Bergdolt 

Prof.  R.  C.  Binder 


DISTRIBUTION  LIST 


No.  of 
Copies 

1 


1 


1 


k 


1 


1 


1 


1 


1 


No.  of 

Copies  Organization 


Organization 


Princeton  University  1 

Princeton,  New  Jersey 
Attn:  Prof.  D.  Ber shader 

Pennsylvania  State 

University  2 

Physics  Department 
State  College,  Pennsylvania 
Attn:  Prof.  R.  G.  Stoner 

University  of  Michigan 
Department  of  Physics  1 

Ann  Arbor,  Michigan 
Attn;  Prof.  Otto  laPorte 

Institute  for  Fluid  Dynamics 
and  Applied  Mathematics  1 
University  of  Maryland 
College  Pai’k,  Maryland 
Attn;  A.  Weinstein 
S.  I.  Pal 
R.  Betchov 
E.  L.  Resler,  Jr. 

University  of  California 
Low  Pressures  Research 
Project 

Berkeley,  California 
Attn:  Prof.  S.  A.  Schaaf 

University  of  Illinois 
Aeronautical  Institute 
Urbana,  Illinois 
Attn:  Prof.  C.  H.  Fletcher 

University  pf  Michigan 
Aeronautical  Research  Center 
Willow  Run  Airport 
Ypsllanti,  Michigan 

United  Aircraft  Corporation 

Research  Department 

East  Hartford,  Connecticut 

Dr.  A.  E.  Puckett 
Hughes  Aircraft  Company 
Florence  Ave.  at  Teal  St. 
Culver  City,  California 


Professor  W.  Bleakney 
Palmer  Physical  Laboratory 
Princeton  University 
Princeton,  New  Jersey 

,  Professor  J.  0.  Hirsch- 
1'  felder 

Chemistry  Department 
University  of  Wisconsin 
Madison  6,  Wisconsin 

f  Professor  H.  W.  Emmons 
Harvard  University 
Cambridge  58,  Massachusetts 

(^’Vprofessor  L.  S.  G. 

Kovaszny 

912  Belgian  Avenue 
Baltimore  18,  Maryland 
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